Magnetic field dependence of the exciton energy in a quantum disk 
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The groundstate energy and binding energy of an exciton, 
confined in a quantum disk, are calculated as a function of an 
external magnetic field. The confinement potential is a hard 
wall of finite height. The diamagnetic shift is investigated 
for magnetic fields up to 40T. Our results are applied to 
In y Al\-yAs/Al x Ga\- x As self-assembled quantum dots and 
very good agreement with experiments is obtained if the light 
hole is assumed to be involved in the exciton formation. Fur- 
thermore, we investigated the influence of the dot size on the 
diamagnetic shift by changing the disk radius. The exciton 
excited states are found as a function of the magnetic field. 
The relative angular momentum is not a good quantum num- 
ber and its value changes with the magnetic field strength. 

71.35.Ji, 71.55.Eq, 73.20.Dx, 31.15.Fx 



I. INTRODUCTION 

Recently, there has been much interest in the study of 
quantum dots, which are structures in which the charge 
carriers arc confined in all three dimensions. Especially 
the self-assembled quantum dots Q are considered to be 
very promising for possible applications, such as quan- 
tum dot lasers due to their large confinement en- 
ergy and high optical quality. The dots are formed by 
the Stranski-Krastanow growth mode in which a mate- 
rial, e.g. In y Al\- y As, is deposited on another mate- 
rial with a substantially different lattice parameter, e.g. 
Al x Ga\- x As ||. The lattice mismatch, which is required 
for this growth process, is typically about 4% Q. Initially, 
the growth is two dimensional, but after a critical thick- 
ness of a few monolayers, coherent islands are formed 
due to strain effects. The shape of the formed islands 
is not well known, but is expected to resemble a lens or 
a pyramid. The density, size and shape of the dots are 
strongly dependent on the growth conditions. Typical 
sizes of dots vary between the basis size of 7 to 20nm 
and a height of a few nanometers. The density of the 
dots is of the order 10 n cm~ 2 M. 

The properties of confined excitons have been the sub- 
ject of many theoretical studies. Bryant || used vari- 
ational and configuration-interaction representations to 
study excitons in quantum boxes. Later, matrix diag- 
onalization techniques were used to study the exciton 
energy in a quantum dot with parabolic confinement po- 
tential. Song et al. || studied the effect of non circu- 
lar symmetric structures, and Halonen et al. |7j] stud- 
ied the influence of a magnetic field. More recently, 



Pereyra et al. || investigated magnetic field and quan- 
tum confinement asymmetry effects on excitons, again for 
the case of parabolic confinement. These studies have 
shown a strong competition between the quantum dot 
size, Coulomb interaction and magnetic confinement. 

In the present work, we approximate the quantum dots 
by a quantum disk with a hard-wall confinement of fi- 
nite height ||, as found in self-assembled quantum dots, 
and include the mass mismatch between the dot and bar- 
rier material. We present a theoretical study of the ef- 
fect of an external magnetic field on the properties of 
an exciton in the quantum disk, fully taking into ac- 
count the Coulomb interaction between the electron and 
the hole. The groundstate energy and binding energy 
of the exciton are studied as a function of the magnetic 
field. This allows us to determine the diamagnetic shift 
of the exciton, which we find in very good agreement 
with the experimentally observed shift by Wang et al. 
]l0| . In most of the previous theoretical work , this dia- 
magnetic shift was only determined for very low values of 
the magnetic field, where the confinement energy is larger 
than the cyclotron energy, and could be approximated by 
e 2 (p 2 ) B 2 /8/j, ]lO|-p^[. In our calculations, we consider 
magnetic fields up to 40T. For such large magnetic fields, 
the weak field approximation is no longer valid, because 
now the cyclotron energy overcomes the confinement en- 
ergy and the particles will act rather as free particles in 
a magnetic field Jl3[] . We find that the magnetic field 
dependence of the diamagnetic shift can be very closely 
approximated by (}B 2 /(l + aB). To be able to make 
a valid comparison between theory and experiment, we 
considered for our simulations Jno. 55^0. AbAs quantum 
dots, which were experimentally studied by Wang et al. 

The paper is organized as follows. In Sec. II, we 
present the theoretical model and explain our method 
of solution. The results for the exciton groundstate and 
the comparison with the experimental results of Ref. E^] 
are presented in Sec. III. In Sec. IV, we describe the 
effect of changing the disk radius on the exciton energy 
and diamagnetic shift. The results for the exciton en- 
ergy spectrum are presented in Sec. V. Our results are 
summarized in Sec. VI. 



II. THEORETICAL MODEL 



The Hamiltonian describing our system is given by 
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where the indices j — 1,2 correspond to the electron and 
the hole with masses mi,m2, respectively, Vj{pj,Zj) — 
0(pj < R,\zj\ < d/2),Vj yQ (otherwise) is the confine- 
ment potential with R the radius of the quantum disk 

and d its thickness, pj — \J + y], V c (r) — — e 2 /e|r| , 

and qj = =pe. Here and below the upper and lower sign 
correspond to electron and hole, respectively. For con- 
venience we will sometimes also use the notations e, h 
instead of 1,2. We allow for a difference in mass be- 
tween the dot region and the region outside the dot: 



m 3 -(r) 



rn„ 



inside the disk and uij(r) — mtj out- 



side the disk. In our numerical work, we used the fol- 
lowing values for the physical parameters: e = 12.71, 
m w<e = 0.076m , m&, e = 0.097to , m Wt h = m bth = 
0.45m , V e , — 258meV, and Vh, a = 172meV, which are 
typical for the Ino.55Ala.45As/Alo.35Gao.65As system. 

Using cylindrical coordinates r,j — (zj, pj, <j)j) the one- 
particle Hamiltonian takes the form 
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dzj m,j dzj pj dpj m,j dpj Pj m j dcp 2 I 

i d 1 

To^JlT + + V i( z j>Pj) > ( 3 ) 



where 1 



eB/rrijC are the electron and hole cyclotron 



frequencies and the vector potential is taken in the sym- 
metrical gauge A =\Bpe ( f ) . 

The one-particle wave functions are separable 
*j(z,p,<l>) = (l/v^e^&fo.p,-), where I = 
0, ±1,±2, ... is the angular momentum, and the wave 
functions £, l ji(zj,pj) are eigenfunctions of the Hamilto- 
nian 
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where the index i denotes the eigenenergies of flj. As a 
consequence of the axial symmetry of our problem, there 
is no coupling between the wave functions with different 
values of the total angular momentum L. Therefore, we 
can construct the exciton wave function with fixed 
total momentum L as the linear combination 



with E the eigenenergy, for brevity x denotes the co- 
ordinates (z%, Z2,pi.p2), Vq i s the matrix element of the 
Coulomb interaction 

e 2tt 
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L m = 2l m + 1 is the total number of angular harmonics 
in the expansion. 

A common technique to solve the eigenvalue problem 
is to use an expansion of the wave function in a suitable 
set of basis functions. For the typical sizes of the quan- 
tum disks considered here, the exciton binding energy is 
much smaller than the confinement energy. As a con- 
sequence, a natural choice is to take the eigenfunctions 
A i of the one-particle Hamiltonian. But for our present 
problem such an approach runs into obstacles because 
of the enormous number of basis functions which are re- 
quired to obtain the binding energy with sufficient ac- 
curacy. Indeed, using the one particle eigenfunctions for 
different values of the angular momentum I and quantum 
number i = 1, one has to calculate L m I A matrix el- 
ements of the Hamiltonian. In the present case of hard 
wall confinement, the one dimensional eigenfunctions are 
too complicated in order to obtain an analytical expres- 
sion for the Coulomb matrix elements. Therefore, a nu- 
merical integration procedure has to be used on the space 
grid with size N g = (K x A^) 2 , where K, N are the num- 
bers of grid points for the longitudinal and transverse 
directions, respectively. In principle, the difficulties in 
the calculation of the Coulomb matrix elements can be 
avoided by applying an appropriate basis, for instance 
the nonorthogonal Gaussian basis, which is widely em- 
ployed in quantum chemical simulations. But in this case 
there is an increase of the number of functions, which are 
needed, leading to difficulties with diagonalizing a large 
L m I 2 x L m I 2 non-sparse matrix. Note that for an arbi- 
trary basis, the number / — i z x i r is determined by the 
number of one-particle wave functions in the longitudinal 
(i z ) and the radial (i r ) directions. The total number of 
operations depends crucially on the considered number 
of subbands i z in the z-direction. For a small ratio d/R 
of the longitudinal to transverse size of the quantum disk 
as given before, we can limit ourself by taking only one 
subband Hl5|. 
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where the functions ip l (x) obey the Schrodinger equation 



A. 3D exciton problem 

For arbitrary values of the ratio d/R we present a nu- 
merical technique based on the use of a finite difference 
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scheme. Let zu, {k — 1, K), p n , (n — 1, N) be some 
nonuniform space grid in the longitudinal and transverse 
directions for both electron and hole coordinates. Using 
the appropriate symmetry conditions for the ground wave 
function in the longitudinal direction dif)/dzj(zj=0) = 
we can limit ourselves to the region Zj > 0. Thus, the 
first point of the z-grid corresponds to z = 0. The up- 
per [zk > d/2) and right {pn > R) boundaries of the 
simulation region correspond to the barrier region where 
the wave function and its derivatives go to zero. There- 
fore, the Neumann conditions dtp/dzj = 0, dip/ dp j = 
are employed for these boundaries. To obtain the finite 
difference scheme for the one-particle Hamiltonian, in- 
cluding the discontinuous behavior of the particles mass 
and external potential, we integrate the expression over 
the square (z k _ 1/2 < z < z k+1/2 , p n -i/2 < P < Pn+1/2), 
where the subgrids with noninteger indexes are deter- 
mined by the relations z k+l/2 = (z k+1 + z k )/2, p n+1 / 2 = 
(Pn+i + Pn)/2, 2-1/2 = P-1/2 = 0. Substituting the fi- 
nite difference expressions for the derivatives of the wave 
function dip/dz{z = z k+1/2 ) = {tpk+i - ipk)/(zk+i - 

Z k ),dlp/dp{p = p n+ i/ 2 ) = {tpn+l - 1pn)/(Pn+l - Pn) we 

obtain the following finite difference scheme for the one- 
particle Hamiltonian 



h z,i Iz" + !^ m i 1 ( z > P=Pn)dz, 
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with the coefficients 
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,k,n=N n 
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J z.k p,nJz k _ 1/2 Jp n - 1/2 Jy r 

Once a 
difference Hamiltonian H = 5i t i> Y^j=i Rj + Vc' 1 $M,M' 
has been constructed, we have to develop a technique to 
obtain the ground state of the sparse matrix R . Here Sij 
is the unit matrix, index M denotes all indexes corre- 
sponding to the space grid. Note that the number of non 
zero elements of the matrix H is only proportional to 
L% n N g . This is a key distinction from the commonly ac- 
cepted expansion over basis functions, where this number 
increases as the second power with the number of func- 
tions. However, the size of our matrix is still large and 
therefore direct diagonalization methods are not suitable 
for solving our problem. The best suitable approach to 
find only the lowest eigenvalue E g and eigenvector ^ is 
the inverse iteration method, where the eigenvector ty l 
at the i th stage of the iteration is obtained by solving the 
following equation 



{H-\8i tl ,8 M , M ,)W = ¥~ 1 
with the subsequent normalization 
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where the brackets ( , ) stand for scalar multiplication. 
The eigenenergy is obtained in the usual way E l g = 

(^ % ,H i $> % \ . The value of the parameter A < E g is cho- 
sen such that a minimum absolute value of the matrix 
(H — \6iji6m,m') corresponds to the ground state of the 
matrix R. There exist many numerical relaxation tech- 
niques to solve the boundary value problem. Using stan- 
dard methods one has to solve the equation with good 
precision at each stage of the inverse iteration procedure. 
Here, we propose a new technique, which generalizes in 
fact the commonly accepted Gauss-Seidel methods |l|] 
with inverse iterations. The value of the eigenvector W 
for the mesh points (I, m = fei, rti, k 2 , n 2 ) is obtained by 
using the following relation 

2 

*i - + ai-iGx + 6 2 )/ C^Pj + U c " - A) , (12) 



where h z> k = z k+1/2 - z k _ 1/2 , h p , n = {p 2 n+l/2 - 
Pn-1/2) ^ ue ^° discontinuity of the mass and 
the external potential at the disk boundary, special 
care must be taken in the choice of the expression 
for its grid values. In the expressions the aver- 

j i r i-i k.n k,n k,n j 

aged value ot tne masses m^ z , nij' , rrij and poten- 
tial Vj ' n are determined by the following relations 

(mf")- 1 = hp.nfp^j* pmj 1 (z^z k ,p)dp, (m^")" 1 = 



where 



-E^^'Ci. (13a) 



(13b) 
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For the ground state 'J'i = ^>i-i/(E g — A), we found that 
the maximum rate of convergency is realized by using 
the following values of the parameters = X/{E % g — A) 
and Q!j = i = 1. 



B. 2D exciton problem 

For quantum disks with large radius i? 3> d we use the 
adiabatic approach, a technique which was already suc- 
cessfully applied in Refs. |l4],|l7|]. Within this approach, 
we can write the wavefunction as 



(14) 



where ijij (Zj ) corresponds to the groundstate of the lon- 
gitudinal Hamiltonian 



H 



z,3 



t?_d_J_J!_ 

' 2 dzimtdzi jAZj) 



(15) 



for electron and hole (j=2), respectively. Since 

the wave function penetrates only slightly into the bar- 
rier region in the radial direction, the longitudinal be- 
havior of the effective masses rrij and the confinement 
potentials Vj tZ can, to high accuracy, be approximated 
by rrij tZ — rrij(z,pj = 0), Vj. z — Vj(z,pj = 0). Then 
the wave function of the ground state has a simple 
form inside, ipj(\z\ < d/2) — cos(kjz), and outside, 
■0j(|z| > d/2) = exp(— Kj \z\), the disk, where kj = 
yf2m^jEo~^j/K and Kj = \/2mb,j(Vo,j — Eo <z j)/h. The 
energy of the groundstate Eo tZ j is obtained from the con- 
tinuity of the wave function and conservation of the cur- 
rent to _1 9'0/9z at the boundary (|z| = d). Substitut- 
ing expression ([l4]) into the Schrodinger equation and 
integrating out the Zj coordinates by taking the average 
(il>i(zi)ip2(z2)\H\ipi(zi)ip2(z2)) = H 2D , we obtain the 
effective two dimensional Hamiltonian 
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obtained an analytical approximation to the effective 
Coulomb potential 



K'(p) 
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where K$ is the modified Bessel function. For a sys- 
tem with infinite barriers the value 7 = 0.277<i gives 
the best fit to the effective Coulomb potential. As 
a consequence of the penetration of the electron and 
hole into the barrier region, the value of 7/d increases 
with decreasing disk thickness. We have found that for 
our parameters of the quantum disk, R = 8.95nm and 
d = 3.22nm, the value 7 = 1.675nm gave the best 
fit to the results obtained from direct numerical cal- 
culation of the effective two-dimensional Coulomb po- 
tential. The other parameters of the two-dimensional 
Hamiltonian are Eo, ze = 116.06meF, -Eo,z/i = 38.13meV, 

133.87mey, m' we = 0.080m , 
= 0.45mo, where indices 



VU = 141.94meV, V' oh 
m' be = 0.097to , m' wh 



'Mi 



+Vc(pl - P2) 



(16) 



where pj = -ihd/dpj, V-{p 3 ) = Vj(zj = 0,pj) - E , z j, 
irij(pj > R) = m bJ , 



e, h correspond to electron and hole, respectively. The 
numerical diagonalization technique for the 2D Hamilto- 
nian was presented already in Ref. E4] . 



III. RESULTS AND DISCUSSION FOR THE 
EXCITON GROUNDSTATE 

We have calculated the exciton groundstate energy and 
exciton binding energy as a function of an applied mag- 
netic field. We used for our simulations the physical pa- 
rameters of the Ino.55Alo.45As self-assembled quantum 
dots, used in the experiment by Wang et al. [^0|. The 
studied disks have a height of 3.22nm and a radius of 
8.95nm. The other parameters were already given above. 
Fig. 1 shows the probability distribution of the electron 
(solid curves) and hole (dashed curves) |^i(ri)| 2 ,i = e,h, 
for the ground state along (z = 0, p) and perpendicular 
(z,p = 0) to the disk under consideration. Along the p- 
direction, the electron and hole are confined within the 
disk but along the z-direction, there is appreciable pene- 
tration into the barrier material. 

In Fig. 2 the exciton groundstate energy is plotted as a 
function of the magnetic field. This groundstate energy 
is given by 



1 



1 



d/2 



dzj \ipj(zj)\ 



m'(pj < R) m w j J Q 

1 f 00 

m b,j Jd/2 

and the effective Coulomb interaction is 

o2 r°° 



E =E e + E n + E e 



(20) 



e z f°° 

Vc(p) — / dzidz 2 

e J-00 



iy*i)l 2 hM* 2 )| 2 

y/( Zl -Z 2 )2 + \p\ 2 



(17) 



(18) 



Using a Gaussian shape for the longitudinal wave func- 
tion of the groundstate, the authors of Ref. Jl8| have 



where E e and E h are the single electron and hole ener- 
gies, respectively and E exc is the exciton binding energy. 
The solid curve shows the result of the full 3D treatment 
of the problem, whereas the dashed and dotted curves 
are calculated using the adiabatic approximation. For 
the latter case, we make a distinction between the cases 
with and without correlation. For the case without cor- 
relation, the Coulomb interaction is calculated using the 
single particle electron and hole wave functions 
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E. 



lfj e ifj h 



(21) 



The total exciton wave function was used in order to 
calculate the energy with correlation 



E, 



e . h 



1 



e . h 



(22) 



Figure 2 shows an enhancement of the groundstate en- 
ergy with increasing magnetic field for all three cases. 
The correlation energy, which is given by the difference 
between the dotted and dashed curve, is 3AmeV for 
B = OT and increases to AAmeV for B = 40T. 

The inset shows the exciton binding energy as a func- 
tion of the magnetic field. Again we see an increase for 
increasing magnetic field as expected. This is not sur- 
prising, because by applying higher magnetic fields the 
particles are more confined, they are closer to each other 
and therefore more tightly bound, which implies an in- 
crease of the binding energy. The assignment of the dif- 
ferent curves is the same as for the groundstate energy 
(see main figure). Note that the inclusion of correlation 
increases the binding energy at B = OT with 14.5% while 
the full 3D treatment of the problem further increases the 
binding energy with 13.6%. 

From our calculation of the exciton groundstate en- 
ergy, we can easily determine the diamagnetic shift of the 
exciton, which is defined by AT = E{B) - E(B = 0). 
The result is shown in Fig. 3, where the curves indicate 
our calculated results for the three cases, as mentioned 
above, and the squares are the experimental results, as 
obtained by Wang et al. From the comparison be- 

tween the different approaches and experiment, we no- 
tice: 1) for B < 8T all three approaches give practically 
the same result which agrees perfectly with experiment, 
2) when B is increased above 8T the three theoretical 
approaches have the same qualitative independence but 
there are small quantitative differences in the slope of the 
curves, and 3) in the high field regime, i.e. B > 20T, 
our theoretical results substantially underestimate the 
experimental result. The masses used for these calcula- 
tions were the ones given by Wang et al. in Rcf. |T(]] 
(m Wl e = 0.076mo, = 0.097mo: m h = 0.45mo) and 
it is clear that here the heavy-hole mass was used. How- 
ever, in Ref. jl9| it was argued that for a magnetic field 
normal to the sample plane, the light hole mass should 
be used. Because the dot height is much smaller than 
the dot radius, heavy hole character is expected in the 
growth direction for the ground hole state and light hole 
character for in-plane motion. Therefore, for B normal to 
the sample plane, the light hole mass should be used. In- 
cluding the effects of strain, they find for InAs dots that 
m e — 0.055mo and rrih — O.lmo- Combining this with 
values for AlAs (20), we find by linear interpolation to 
the material Ihq ^AIq ^As values of 0.080mo and 0.2mo 
for respectively the electron and the hole mass. The re- 
sult for the diamagnetic shift in this case is depicted as 



the dot-dashed curve in Fig. 3. This result is in very good 
agreement with the experimental results. 

In previous theoretical work, only the exciton energy 
and wave function at B = were considered, from which 
the diamagnetic shift can be calculated as AE = /3B 2 , 
where (3 = e 2 (p 2 ) /8/x and (p 2 ) is the mean quadratic 
electron-hole distance. This is a good approximation 
in case of low magnetic fields, when the magnetic con- 
finement is much lower than the confinement due to the 
quantum dot. However, for higher magnetic fields, the 
magnetic confinement becomes more important. Then 
the quadratic dependence of the energy shift on the mag- 
netic field will change into a linear dependence, due to 
the formation of Landau levels. In this case the energy 
shift becomes AE = h(Lo c>e + u) c ,h)i where oj c .i = eB/rrii 
is the cyclotron frequency. With this knowledge, one can 
construct the function 



AE = 



(3B 2 
1 + aB 



(23) 



which interpolates between the small and large magnetic 
field behaviour and where (3 and a are taken as fitting 
parameters. This formula gives for low magnetic fields 
(B — ► 0) the already known expression AE = (3B 2 , and 
for high magnetic fields (B — > oo) AE = {(3/a)B. It 
turns out that Eq. (23) gives an extremely good fit to 
the numerical results of Fig. 3 for (3 — 6.63p,eVT~ 2 and 
a = 3.25 x lO- 3 T-\ We found that the fitted curve 
reproduces the solid curve in Fig. 3 so well that they 
can not be discriminated. We also calculated (3 us- 
ing the expression (3 = e 2 (p 2 ) /8(jl for B -> 0, which 
resulted into the value (3 = 9.58fieVT~ 2 . This value 
is substantially higher than the one found by fitting. 
In the other limit, we compare AE — ((3/a)B with 
AE = huj c = (he/fi)B, where p = m e m,h/{m e + mh) 
is the effective exciton mass in In AlAs. Such a calcu- 
lation gives he/ fj, = 1.68 x lO^eVT' 1 , which is smaller 
than the fitted value 0/a = 2.04 x lQ^eVT -1 . The fit- 
ted results within the adiabatic approximation with and 
without correlation are respectively, (3 = 7.56fieVT~ 2 , 
0/a = 1.52 x lO^eVT- 1 and = 7.79fieVT~ 2 , /3/a = 
2.75 x 10~ 3 eVT~ 1 . Using the light hole mass instead 
of the heavy hole mass, we find respectively the fol- 
lowing fitted and calculated results: (3 = 9.16[ieVT~ 2 , 
[3/a = 2.27meVT- 1 and [3 = UM^eVT^ 2 , [3/a = 
IMmeVT- 1 . 

In the above calculations we investigated the adiabatic 
shift, which is a relative quantity, and therefore in the cal- 
culation of the groundstate energy, the bandgap was not 
included. But when we want to compare the experimen- 
tal excitation energy, the bandgap of the disk material 
is needed. For B = OT, using the heavy hole mass, we 
found a groundstate energy of E = E e + E h + E exc = 
152.2meV, using Eq. (|2fj|). For the case of the light hole 
mass, which gave a better agreement with the experi- 
mental results, the groundstate energy at B = OT is 
178.5meV. To obtain the total excitation energy, as mea- 
sured in photolumincsccncc experiments, e.g. by Wang 
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et al. , the bandgap energy E g has to be added to this 
equation: 

E = E e +E h + E exc + E g . (24) 

For 

our study, we considered Ino ^AlQ ^As/Alo ^Gao.65-^s 
quantum dots. Without strain, the bandgap energy of 
the dot material was obtained by linear interpolation be- 
tween the result for InAs (E g = 0.41eF) and AlAs (E g = 
3.13eV) which results into E g = l.Q'iAeV, whereas we 
found for the barrier material that E g = 2.083ey @,|§. 
The difference in bandgap between the two materials is 
AE g = 450meV. For the total exciton energy, we now find 
E = 1.634el/+0.152eV(0.1785eVO = 1.786ey(1.8125eV) 
using respectively the heavy (light) hole mass. From Ref. 
fiof , we know that the bandgap difference between the 
dot and the barrier material, corrected for strain effects, 
is AE g = 430meV. This means that the bandgap of the 
dot material has increased with 20meV. For the total ex- 
citon energy, this gives us the final result of E = 1.81eV 
using the heavy hole mass and E = 1.83eV using the 
light hole mass. In the experiments, for B = OT, the 
value of E = 1.894eV was found, which gives a reason- 
able agreement with our theoretical result in view of the 
fact that the composition of the alloy in the dot can, for 
example, not be uniform, the dot size is not known with 
high accuracy, etc. 

Next, we investigated the effect of an applied magnetic 
field on the exciton characteristics, using the parameters 
corresponding to the solid curve in Fig. 3. First we con- 
sidered the one-particle characteristics (z^) 1 ^ 2 , (z^) 1 ^ 2 , 

(Pe) 1 ^ 2 an d (Ph) 1 ^ 2 > where z e , Zh and p e , ph are the 
electron and hole coordinates along the z-axis and in the 
plane, respectively. The results are shown in Fig. 4 and 
were calculated using the full 3D approach. The figure 
shows clearly the squeezing of the exciton due to the 

magnetic field, especially for the in-plane direction. The 

l ji 

mean quadratic electron- hole separations (p^) and 

1/2 

( z eh) §i ve an idea °f the- size of the exciton. We de- 
fined p e h = |~p*e — 7?h\ an d z e h — \z e — Zh \ ■ Notice that 
the size of the exciton is comparable to the disk size. 
We see a more substantial decrease with increasing mag- 
netic field than for the single particle wavefunction, which 
agrees with the increased binding of the exciton. 

In Fig. 5 the percentage of the electron (right scale 
in Fig. 5) and hole (left scale in Fig. 5) wavefunction in 
the dot is shown with varying magnetic field. Both the 
results for the 2D case with correlation and the full 3D 
treatment were calculated. More than 90% of the hole is 
inside the dot while only 71-73% of the electron is inside 
the dot. With increasing magnetic field both the electron 
and hole become more confined inside the dot, indicating 
further the squeezing due to the magnetic field. For the 
hole, we observe a flattening of the curve at very high 
magnetic fields, both for the 3D calculation and for the 
2D case. For high magnetic fields, the hole wavefunction 



is in the p-direction totally confined in the dot. However, 
there is still some extent of the wave function outside the 
dot in the z-direction. But since the magnetic field has 
almost no influence on the z-direction, applying higher 
magnetic fields will not attribute to a further increase of 
the amount of the wavefunction inside the dot and there 
will always be a small part of the wavefunction outside 
the dot. 

Figures 6(a,b) are contourplots of the density distribu- 
tion, of the electron and hole, respectively, along a cross 
section in the middle of the quantum dot perpendicular 
to the y-direction. The electron density is defined as 

\ipe(p e ,z e )\ 2 = J dz h J dp h \ip e ' h (p e ,z e ,p h ,z h )\ 2 , (25) 

and similarly for the hole. The solid curves show the re- 
sult for the case of B = OT, whereas the dashed curves are 
the result for B = 40T. The dashed square indicates the 
position of the disk, which is only one fourth of the actual 
disk size. Due to the magnetic field, we see an increase 
in the density inside the dot, both for electron and hole. 
Along the p-direction the particles become more centered 
in the middle of the dot due to the squeezing by the ap- 
plied magnetic field. However in the z-direction, it seems 
at first sight that there is an expansion instead of the ex- 
pected squeezing, but a closer look (by normalizing the 
function to its central value) tells us that this is not the 
case. Of course, the magnetic field is applied along the z- 
axis and has no direct influence on the exciton behaviour 
in the z-direction. In the p-direction however, the mag- 
netic field brings the electron and hole closer together. 
This implies a stronger interaction and we expect that 
this effect should also be seen in the z-direction. This 
is also the case for the mean quadratic hole (z%) and 

1 /2 

electron-hole separation (z^ h ) (see Fig. 4), which de- 
crease as a function of the magnetic field. 

IV. EFFECT OF CHANGING THE DISK RADIUS 

We investigated the effect of the size of the disk on 
the exciton binding energy which is depicted in Fig. 7 
for the case with the full 3D treatment, for the param- 
eters corresponding to the solid curve in Fig. 3. When 
varying the disk radius from R = lnm up to R = 15 
nm (the dot thickness was fixed to d = 3.22nm), we see 
initially a strong increase of the exciton energy by more 
than a factor 2 and beyond R ~ 2.bnm it decreases slowly 
for increasing R. In the 'large' i?-region, the binding en- 
ergy increases for decreasing disk radius due to the larger 
confinement of the electron and hole wavefunction. The 
electron and hole are forced to sit closer to each other, 
which leads to an enhancement of the binding energy. 
This behaviour continues until the disk radius reaches a 
value of R ~ 2.5nm, where the binding energy reaches 
a maximum value of E exc = A7meV. The decrease in 
the binding energy with decreasing R is due to the fact 



G 



that the wavefunction of the particles start to spill over 
into the barrier material, i.e. the electron and the hole 
become less confined, which leads to a much smaller in- 
teraction and therefore a lower binding energy. This is 
due to the competition between the confinement kinetic 
energy and the barrier material potential energy. This 
is confirmed by Fig. 8, where the percentage of electron 
and hole inside the dot is shown as function of the disk 
radius. For R > 6nm these percentages increase very 
slowly with increasing R. We never reach 100% because 
of the substantial penetration of the wavefunction in the 
barrier material along the z-direction (the thickness of 
the dot is only d = 3.22nm). Note that for R = lnm 
only 1.88% of the electron wavefunction is inside the dot 
but 24.20% of the hole wavefunction. 

The effect of a magnetic field on the dot size depen- 
dence of the exciton energy is also shown in Fig. 7 for the 
case of B = 40T. Notice that the largest B-dependence is 
found for very small and very large R. In both situations 
the confinement of the electron and hole are smallest and 
consequently the ratio between the magnetic energy and 
the confinement energy is largest. For intermediate dot 
size, i.e. 3nm < R < lnm we observe the smallest ef- 
fect of a magnetic field on the exciton energy. This is 
the region of dot size where the confinement potential is 
able to strongly confine the electron and hole to a small 
region in space. 

We also investigated the effect of varying R on the 
electron-hole separation, both in the p and in the z- 

1/2 

direction. Fig. 9 shows the result for (z^ h ) and we 
see a rather high starting value at R = lnm, decreasing 
strongly for increasing R. This high value at small R fol- 
lows from the fact that a large part of the wavefunctions 
is outside the dot, so the particles are not really confined 
anymore, which means that they are farther away from 
each other. The inset of Fig. 9 shows the electron-hole 
separation in the p-direction. Also here, we start with 
a high value at very small R, followed by a strong de- 
crease and a minimum of yj {p 2 eh ) = Z.Ylnm at R ~ 3nm. 
Further increasing R, we find again an enhancement of 

1/2 

(Peh) i which initially is linear in R, but for R > 12nm 
starts to level off and reaches a constant value in the limit 
R — > oo. 

The low magnetic field diamagnetic coefficient (3 — 
e 2 (p^h) /(8/u), can be obtained from the results of Figs. 8 
and 9, where fi is the effective exciton mass. The result 
is shown in Fig. 10 and we see a similar behaviour as 
for the radial electron- hole separation (inset of Fig. 9). 
When calculating (3, we took into account the variation 
of the effective exciton mass fi with varying disk radius. 
The effective mass is defined as 

- = — + —, (26a) 
fx m e m h 

with 



m e m e , b m e , w ' 



where m e ^ = 0.097mo and m e w = 0.076mo are the ef- 
fective electron masses in respectively the barrier and the 
well and P w is the probability to find the electron in the 
well. In Fig. 8 we showed that there is a considerable 
change of P w for varying R, and this will have an effect 
on m e and fi. For the hole we have the same mass in 
and outside the well, and therefore m/, = 0.45mo is in- 
dependent of R. In Fig. 11, the evolution of m e and a is 
depicted as a function of R. We see that for very small 
disks, where most of the wavefunction is outside the dot, 
the value of m e converges to m e j, as expected. For larger 
disk radii, this value decreases and for R — > oo, it reaches 
the limit m e = 0.0809mo, which is larger than m ejW due 
to the penetration of the electron along the z-direction 
in the barrier because of the small thickness of the disk. 



V. EXCITON ENERGY SPECTRUM 

The higher radial excited states (N ^ 0), for angu- 
lar momentum L = are calculated within the adia- 
batic approximation. The result for a disk with radius 
R = 8.95nm and thickness d = 3.22nm is shown in 
Fig. 12(a), which clearly shows the appearance of anti- 
crossing of levels for higher N states and the energy scale 
for such states is also substantially larger than for the an- 
gular momentum states. This anti-crossing is due to the 
fact that we consider a fixed angular momentum L for 
all states, which is a conserved quantity. The states with 
fixed L are non degenerate. Again we considered cases 
of different disk radii and we observe an enhancement of 
the anti-crossing for smaller disks (Fig. 12(b)) and a di- 
minishing of the anti-crossing for larger disks (Figs. 12(b) 
and 12(c)). 

To study the anti-crossing more closely, we considered 
the disk with R = 8.95nm and d — 3.22nm (parameters 
corresponding to the dashed/dotted curve in Fig. 3 were 
used) . For the 2D problem, the radial part of the exciton 
wave function for a fixed L can be written as (see also 
Ref. @) 

k=k n n=k n i=im 
TpiPe, Ph) = C l kn R k ^ L+l)/2 (Pe) (27) 

fc=l 71=1 l = -l m 

where k and n correspond to the energy levels of the one 
particle problem of electron and hole, respectively and I 
is the relative angular momentum. The sum ^2 indicates 
that only even values of the relative angular momentum 
I are taken when L is even, and odd values otherwise. By 
studying the values of the coefficients C l kn , we could dis- 
tinguish which one-particle states contribute most to the 
total exciton state. In Figure 13(a), the symbols on the 
curves indicate which is the dominant term, contributing 
to Eq. (27). The inset in Fig. 13(b) gives the (l,n) value 
corresponding to the different symbols. Notice that k re- 
mains 1, while n can have higher values, which implies 
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that the hole excited states are mixed in into the exci- 
ton wave function. The single particle hole states have 
lower energy, due to its higher effective mass. Notice 
that when we connect each symbol by a line, we obtain 
intersecting levels. Such a spectrum would be obtained 
if (I, k, n) would be conserved quantities. Because of the 
electron- hole interaction the different (l,k,n) single par- 
ticle states are mixed which leads to the anti-crossing of 
the levels. Because L is a conserved quantity, crossings 
between different N states are prohibited and therefore, 
for a fixed N, the system is forced to go to a different 
(I, k, n)-state. 

The symbols in Fig. f3(a) indicate only the dominant 
term in Eq. (27), while the total summation considers 
typically about 750 terms. In Fig. 13(b) we show how 
large the contribution of the dominant term is relative to 
the total sum of all terms. This percentage is defined as 

\C l 1 2 

percentage = 1 fat| — ^ x 10 °- ( 28 ) 

2 \ C kn\ 

k,n,l 

In Fig. 13(b) we only show the result for the case of 
N = 1, 2, 5 in order not to overload the figure. We want 
to emphasize that, as a function of the magnetic field, 
the contribution of the dominant term, which can dif- 
fer with increasing magnetic field, is shown and not the 
evolution of the contribution of a particular state. For 
N = 1, the (0, 1, l)-state appears to be very stable, as 
it stays between 85% and 95%. This means that there 
is very little mixing with other states. The (0, 1, 2)-state 
for N — 2 is also very stable at low fields, but from 
B = 20T, the percentage drops, which indicates that 
another state is becoming important and serious mixing 
occurs. Finally at B = 35T the (— 1, 1, 1)-state becomes 
most important, which can also be seen in Fig. 13(a). 
Now the percentage of the contribution of this (—1, 1, 1 )- 
state is plotted and we see an increase with B. Also 
for N — 5, the transitions between the successive states 
are clearly visible from Fig. 13(b). Each dip corresponds 
to a transition to another state with a consecutive anti- 
crossing of the energy levels. A dip indicates strong mix- 
ing between 2 or even 3 states where the dominant state 
gives only a slightly higher contribution than the other 
important state(s). We see strong dips at B = 12T 
and B — 20T, which indicate the mixing between re- 
spectively (0, 1, 3) < — ► (1, 1, 1) and (1, 1, 1) < — > (0, 1, 3) . 
Comparison with Fig. 13(a) shows that these magnetic 
field values mark also the anti-crossings between respec- 
tively N = 4 < — ► N = 5 and N = 5 < — > N = 6. The 
height of the peaks is an indication of the stability of 
the state. At B = 15T e.g., we see a very strong peak, 
whereas at B = 24T, only a small peak appears. Compar- 
ing with Fig. 13(a) learns that, at the region between 20 
and 30T, there is a strong anti-crossing and, although, 
the system passes through the (0,1,3) state, this state 
never becomes really important. 

In the above discussion, we indicated the states with 



(I, k, n) , where / is an integer number. However, the rel- 
ative angular momentum I is not a good quantum num- 
ber, and therefore the expectation value of the operator 
l z = (h/i) d/d — fa) is expected not to be an inte- 
ger. In Fig. 14 the expectation value (l z ) is depicted for 
the different A-states as a function of the magnetic field, 
where (l z ) is calculated by 

<u = E*l c U 2 - ( 29 ) 

l,k,n 

Notice that (l z ) tends to approach an integer value / when 
one of the terms in the sum of Eq. (27) dominates. The 
transition between states with different I is continuous. 
The more stable a state is, the better it approaches an 
integer value I. The result for A = 7 e.g. starts from 
(l z ) = at B = 0T, then increases up to about 0.7571 
and at B = 10T, drops down to -1.8S until B = 20T 
where it starts to increase again, more slowly now, up 
to (l z ) = 0.9% for B = 25T until finally at B = 40T it 
drops to less than —2.5%. This agrees very well with the 
predicted integer values for I in Fig. 13(a). For other N- 
states, the agreement might be less good, which is due to 
the higher mixing with other states. 

Finally, we considered the energy states for different 
values of the total angular momentum L within the adi- 
abatic approximation. The result for a disk with ra- 
dius R — 8.95nm and thickness d = 3.22nm is shown 
in Fig. 15(a) for the lowest radial state N = 1. Notice 
that for B = 0T, the states with L and — L are degener- 
ate, which is lifted by a magnetic field. The correspond- 
ing splitting is the well-known Zeeman splitting. For a 
smaller disk radius, R = 5nm, all energies are shifted to 
higher values, the splitting between the energy levels is 
larger, and the Zeeman splitting is increased (Fig. 15(b)). 
When increasing the disk radius R, i.e. R = 15nm and 
R = 30ram, the difference between the different angu- 
lar momentum levels decreases and the energy shifts to 
lower values (Figs. 15(c,d)). As in previous case of differ- 
ent A-states, also here we studied which single particle 
(I, k, n)-states are most important in the sum of Eq. (27) 
and the percentage of their contribution. Fig. 16(a) de- 
notes, for a disk with radius R = 8.95nm and thickness 
d = 3.22nm, the energies of the different L-states and the 
symbols indicate which (/, n)-state is most important at 
a particular value of the magnetic field. Note that here 
both k and n remain 1. In Fig. 16(b) the percentage 
contribution of the particular (l,k,n) state is depicted. 
We see a transition occurring for the L = —1,-2 and 
—3 states. This follows also from Fig. 17, where the ex- 
pectation value (l z ) of the relative angular momentum 
operator is plotted. For the L = 0,1,2,3 states, (l z ) 
remains quite constant, whereas for the other i-states, 
(l z ) decreases towards a lower value of I. Because the to- 
tal angular momentum is a conserved quantity, energy 
levels corresponding to different L- values are allowed to 
cross, they do not mix. 
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VI. CONCLUSIONS 

We calculated the groundstate energy (and the excited 
states), the binding energy and the diamagnetic shift of 
an exciton in a quantum disk with radius R and thickness 
d for a hard wall confinement potential of finite height. 
The mass mismatch between the dot material and the 
surrounding material was taken into account. Our calcu- 
lation is based on the finite difference technique, where we 
used three different theoretical approaches, which include 
the electron-hole correlation on different levels. The 3D 
treatment is valid for arbitrary values of R and d and 
provides an 'exact numerical' treatment of the exciton 
problem. For R 3> d, the adiabatic approach is applica- 
ble and here we distinguish the cases with and without 
correlation. The latter only uses the single particle wave 
functions in order to calculate the exciton binding energy, 
whereas the first uses the total exciton wavefunction. 

Under the influence of an external applied magnetic 
field up to 40T, we find an increase of the exciton ground- 
state energy and binding energy. The electron-hole sep- 
aration shows a squeezing of the exciton due to the mag- 
netic field. This can also be seen from the electron and 
hole densities in and around the dot. Our theoretical re- 
sults of the diamagnetic shift are in ver y go od agreement 
with the experimental results of Ref. |f0[ if we assume 
that the light hole is involved in the exciton. 

When considering a varying disk radius R, we found 
a strongly decreasing exciton binding energy with de- 
creasing R for very small i?-values, which indicates that 
the dots are too small to confine the exciton. This ex- 
planation is corroborated by an investigation of the ra- 
dial electron-hole separation and of the percentage of the 
wavefunction in the dot, which indeed shows that, for 
very small R, a large part of the wavefunction is situated 
outside the dot. In the large R- regime the exciton bind- 
ing energy decreases with increasing R and approaches a 
constant value for R — > oo. In the presence of an applied 
magnetic field, the exciton binding energy approaches a 
constant value for large disks much earlier than for the 
B = OT case, indicating that the dot confinement is dom- 
inated by the magnetic confinement. 

Results for higher excited radial states, N > 0, show 
an anti-crossing of levels which is more pronounced for 
small dot radius. The total angular momentum L is a 
conserved quantity. The relative angular momentum I, 
however, is not a good quantum number. Because of the 
coupling between the electron and the hole, the exciton 
wave function is a linear combination of all possible one- 
particle wave functions. We investigated which (I, k, re- 
states contribute most and how large its contribution is 
to the total exciton wavefunction. Furthermore we in- 
vestigated the expectation value of the relative angular 
momentum operator l z , which is not quantized and varies 
with the magnetic field. The degeneracy of the different 
total angular momentum states is lifted due to the pres- 
ence of the confinement potential and the Zecman split- 



ting. This splitting decreases with increasing dot radius 
R. Also here we investigated the contribution of the one 
particle states to the total exciton wavefunction. The en- 
ergy states with different total angular momentum L can 
cross with varying magnetic field, because L is a good 
quantum number. 
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FIG. 1. Side view of the quantum disk together with the 
electron (solid curves) and hole (dashed curves) probability 
distribution along the (p, z = 0) and the (p = 0, z) direction. 

FIG. 2. The exciton groundstate energy as a function of 
the magnetic field. The solid curve shows the result obtained 
within the full 3D treatment, whereas the dashed and dotted 
curves are the result for the adiabatic approximation, respec- 
tively with and without correlation. In the inset, the exciton 
binding energy is plotted. The same curve conventions are 
used as in the main figure. 

FIG. 3. The diamagnetic shift of the exciton energy as a 
function of an external magnetic field. The curves are our 
theoretical results within different approximations and the 
squares are the experimental results of Wang et al. fllfj] . 



FIG. 9. Electron-hole separation in the z-direction as a 
function of R. The inset shows the radial electron-hole sepa- 
ration as a function of R. 

FIG. 10. The diamagnetic coefficient f3 as a function of the 
disk radius R for d — 3.22nm. 

FIG. 11. The effective electron (solid curve) and exciton 
(dashed curve) masses as a function of the disk radius R. 

FIG. 12. The exciton energy states for different values of 
the radial quantum number N, as a function of the mag- 
netic field for L = 0. Four different disk radii are considered: 
R — 8.95nm (a), R = 5nm (b), R = 15nm (c) and R — 30nm 
(d). 

FIG. 13. (a) The higher excited radial states for 
R — 8.95nm and d = 3.22nm. The symbols indicate the (Z, n) 
value of the most important single-particle state, where I is 
the relative angular momentum and n is the radial hole quan- 
tum number, (b) The percentage of the contribution of the 
most important state as a function of the magnetic field. 

FIG. 14. The expectation value of the relative angular mo- 
mentum operator Z z as a function of the magnetic field, for 
L = and different iV-states. 



FIG. 4. The extent of the electron, the hole and the exciton 
as a function of the magnetic field, along the radial (a) and 
longitudinal (b) direction, respectively. 

FIG. 5. The percentage of the wavefunction in the dot as 
a function of the magnetic field, both for the electron (right 
axis) and the hole (left axis). The solid curves are the result 
for the 2D case with correlation, the dashed curves are the 
result obtained within the full 3D treatment. 

FIG. 6. Contourplot of the (a) electron and (b) hole den- 
sity in a plane through the center of the quantum dot with 
size R — 8.95nm and d — 3.22nm. The plot shows only one 
quarter of the total space. Results are shown for B — 0T 
(solid curves) and for B = 40T (dashed curves). 



FIG. 15. The exciton energy states for different values of 
the total angular momentum L, as a function of the magnetic 
field. Four different disk radii are considered: R = 8.95nm 
(a), R = 5nm (b), R = 15nm (c) and R = 30nm (d). 

FIG. 16. (a) The energy states for different total angular 
momentum L for a disk with R = 8.95nm and d — 2>.22nm. 
The symbols indicate the dominant single-particle (Z, n) states 
and (b) gives the percentage of the contribution of this term 
to the total wave function. 

FIG. 17. The expectation value (Z z ) for the different 
Instates as a function of the magnetic fields. 



FIG. 7. The exciton binding energy as a function of the 
disk radius R for a disk thickness of d = 3.22nm. Results are 
shown for two magnetic fields as indicated. 



FIG. 8. The percentage of the electron (solid curve) and 
hole (dashed curve) wavefunction in the dot as a function of 
the disk radius R. Symbols are the calculated points and the 
curve is a guide to the eye. 
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